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^ \ The polarized thermal radiation at far infrared and submillimeter wavelengths 

\ from dust grains in interstellar clouds with irregular magnetic fields is simulated. 

^ \ The goal is to determine how much irregularity in the magnetic fields can be 

\C) \ consistent with the observations that the maps of the polarization vectors are 

^ \ relatively ordered. Detailed calculations are performed for the reduction in the 

\ fractional polarization and for the dispersion in position angles as a function of 

Tjj- \ the ratio of the irregular to the uniform magnetic field and as a function of the 

^ \ relevant dimensions measured in correlation lengths of the field. We show that 

'-^ ' the polarization properties of quiescent clouds and of star-forming regions are 

^ \ consistent with Kolmogorov-like turbulent magnetic fields that are comparable 

\ in magnitude to the uniform component of the magnetic fields. If the beam size 

^ \ is much smaller than the correlation length Lcorr of the fields, the calculated 

\ percentage polarization p decreases to an asymptotic value when the number of 

\ correlation lengths Neon through the cloud exceeds a few tens. For these values 

\ of NcoTT^ the dispersion in the position angles a a is still appreciable — decreasing 

■ - - ' to only about 20°. However, when the finite size of a telescope beam is taken into 

account, the asymptotic value of p is reached for fewer correlation lengths (smaller 
Neon) due to averaging over the beam; a a becames much smaller and consistent 
with the observational data. The smoothing of the polarization properties due 
to the combined effect of the thickness of the cloud and the finite size of the 
beam can be described by a single variable which we designate as the generalized 
number of correlation lengths. 

In addition, we study various factors that may contribute to the decrease 
in the linear polarization percentage with increasing intensity that is observed 
at submillimeter and far infrared wavelengths in many, though not in all, dark 
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clouds (the "polarization hole" effect). Depolarization due to a density cutoff 
in the polarizing effect of dust, due to thermalization, and due to correlations 
between the density and the properties of the magnetic field are considered. 

Subject headings: dust, extinction — ISM: clouds — ISM: magnetic fields — 
MHD — polarization — turbulence 

1. Introduction 

Since its discovery in late 1940's (e. g., Hiltner 1956), the polarization of starlight 
through selective absorption by interstellar dust grains has been considered as a valuable 
tool to study both the large-scale and small-scale structure of interstellar magnetic fields. 
However, about ten years ago it was realized that the polarization of starlight may not be a 
good probe of the magnetic fields in dense interstellar clouds. 

The polarization is caused by aligned grains and it is natural to expect that the frac- 
tional polarization increases when the light from a background star passes through a region 
of higher extinction. However, the observations give a different picture. Neither the per- 
centage polarization, nor the pattern of the polarization seem to differ appreciably for rays 
passing through dense clouds in comparison with the rays that pass through the neighboring 
intercloud medium (Myers & Goodman 1991). Arce et al. (1998) demonstrated that the 
empirical relationship between the maximum percentage polarization of starlight Pmax and 
the optical extinction Ay 

Pmax(%) = 3^v (1) 

is only valid for stars that lie behind the "general" interstellar medium and does not hold 
for stars that arc located behind dense clouds. They concluded that in a dense medium for 
which Ay > 1 — 2, the polarizing efficiency of dust grains is significantly reduced due to poor 
alignment or due to changing grain properties (e.g., spherization). That is, the polarization 
of starlight cannot be a probe of the magnetic field structure in dense quiescent and star- 
forming regions. This would also explain why there is little correlation between the starlight 
polarization pattern and quiescent cloud structure (Myers & Goodman 1991) which would 
be expected in the medium dynamically controlled by the magnetic field. 

However, the observed inefficiency of the polarization of starlight in dense clouds is at 
variance with the growing body of observations of the polarized thermal emission by dust 
(e.g. Hildebrand et al. 2000). Significantly, this emission is observed not only in regions of 
ongoing star formation where the conditions for alignment can be different from those in dark 
clouds, but also in dark and apparently starless cores (Ward-Thompson et al. 2000; Crutcher 
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et al. 2004; Jones 2003). According to Ward-Thompson et al. (2000), the polarization of 
the emission at 850 /im by cold (~ 10 K) dust reaches 10% at the L1544 and L183 cores. 
These cores have no internal sources of radiation, and their external illumination should be 
negligible as some molecules (that arc easily destroyed by UV radiation) are present even in 
the outermost parts of these cores, well outside of regions where the polarized emission is 
observed (Dickens et al. 2000; Tafalla et al. 2002). This implies that Ay > 3 even at the core 
peripheries; the extinction is greater in the inner parts of the cores {Ay > 3). Hence, dust 
apparently preserves its polarizing properties at much higher mean extinctions then those 
probed by starlight polarization. On the other hand, there is recent evidence that a very 
low-mass protostar is present in the L1014 core that previously was believed to be starless 
(Young et al. 2004). This result implies that faint radiation sources may be embedded in 
some other "starless" cores as well, though such sources are presently undetected. 

In our previous paper (Wiebe & Watson 2001) we performed calculations for the hnear 
polarization of starlight due to extinction by aligned dust grains when the starlight traverses 
a medium with irregular magnetic fields. We found that the polarization properties of the 
starlight — the average fractional polarization and the dispersion in position angles — can be 
essentially unchanged if the rms of the irregular component in the optically thick medium 
is greater than the average magnetic field. Thus, the observed tendency for the starlight 
polarization to remain nearly constant as the light passes through a dense cloud is understood 
as due to the changing properties of the magnetic field without the need for a decrease in 
the polarizing efficiency of dust. In this paper we explore the implications of this suggestion 
for the polarized emission by dust, and in particular, how much irregularity in the magnetic 
fields is consistent with the relatively ordered patterns that are observed for the polarization 
directions of the emitted radiation. 

The breadths of spectral lines from interstellar clouds show velocity dispersions that 
exceed the thermal velocities. This and other evidence indicates that turbulent or wave 
motions are pervasive in the interstellar gas (e.g., Elmegreen & Scalo 2004 for a current 
review). From MHD considerations, it follows that irregularities in the magnetic fields are 
expected to be associated with such motions of the gas (e.g., Vazquez-Semadeni et al. 2000). 

The possible origin of polarization holes is also investigated. In general, a polarization 
hole represents a factor of a few decrease in percentage polarization as the intensity increases 
by an order of magnitude, observed in many, though not in all, dark clouds. Its name stems 
from the fact that it is mainly observed in centrally peaked sources, though the anticorrelation 
of the percentage polarization and the intensity may be independent on the source geometry 
(Matthews et al. 2001). This effect is wide-spread but not ubiquitous. For example, it 
does not show up in the NGC 7538 region observed by Momose et al. (2001) at 850 iim. 



-4- 



Akeson et al. (1996) observed the young stellar object NGC 1333/IRAS 4A with high 
angular resolution (5") and found polarized (4%) dust emission arising from very dense gas 
(n > 10* cm~^). An even higher percentage of polarization (10%) is found in the very center 
of the CB68 globule by Valee, Bastien, & Greaves (2000). The conspicuous example of the 
ambiguity about the polarization hole is represented by the Klcinmann-Low Nebula in the 
OMC-1, where the hole is present when observed with low angular resolution at 100 iim 
(Schleuning 1998) and disappears at higher resolution and larger wavelength (Rao et al. 
1998). All this diversity seems to suggest that the polarization hole effect is not related to a 
single mechanism but instead results from a combination of several mechanisms. In addition 
to the poor alignment of grains and the loss of their polarizing ability, other contributing 
factors may include specific variations of the large-scale magnetic field and its unresolved 
small-scale structure that may become more complicated in dense (possibly collapsing) parts 
of the star-forming region. 

We first examine in detail how the irregularities in the magnetic field within a cloud 
reduce the fractional polarization and cause dispersion in the position angles of this polariza- 
tion. That is, we calculate the "polarization reduction factor" F and the rms of the position 
angles function of the ratio of the irregular component to the uniform component of 

the magnetic field, the dimension of the cloud along the line-of-sight A'corr measured in cor- 
relation lengths of the field, and the beamsize of a telescope measured in correlation lengths. 
The only expression in the literature for F of which we are aware (Lee & Draine 1985) 
is applicable only in the limit for asymptotically large A^con and negligible beamsize (even 
with these restrictions, we find that it is inaccurate for certain ratios of the irregular to the 
uniform component of the magnetic field). An expression is available for function of 

Neon (Myers & Goodman 1991), but only for negligible beamsize. Heitsch et al. (2001) have 
calculated the dispersion in polarization angles of dust emission caused by irregular magnetic 
fields for specific, turbulent MHD models. Although beamsize effects were included there, 
the focus of that study (the Chandrasekhar-Fermi relationship) was somewhat different from 
our investigation. 

The paper is organized as follows. In § 2 we describe the methods that are used to 
generate the turbulent magnetic fields as well as the methods to calculate the emission of the 
polarized radiation from the cloud. In § 3, the polarization characteristics of the emergent 
radiation are calculated as a function of the relevant parameters for the case of constant 
density. Effects of the non-uniform density are outlined in § 4 where the emphasis is on 
considerations that might contribute to the polarization hole effect. A concluding discussion 
is given in § 5. 
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2. Basic Methods 

We describe the magnetic field as a sum of uniform and turbulent components. As 
in our previous investigations, representative turbulent magnetic fields are created by sta- 
tistical sampling of the Fourier components of a power spectrum with a Kolmogorov-like 
(i.e., power law) form and with Gaussian distributions for the amplitudes (Wallin, Watson, 
& Wyld 1998, 1999; Watson, Wiebe, & Crutcher 2001). Such methods are standard (e.g., 
Dubinski, Narayan, & Phillips 1995). As in our previous investigations, we focus on a power 
spectrum that is somewhat steeper than Kolmogorov (here, wavenumber k~^f^ instead of 
the Kolmogorov fc~^/"^) based in part on the values of the correlation lengths that we extract 
(see Watson et al. 2001) from the results of MHD simulations by others (Stone, Ostriker, 
& Gammie 1998). Two quantities characterize the turbulent magnetic fields created in this 
way — the rms value of each of the spatial components of the irregular magnetic field -Brms 
(assumed isotropic) and the correlation length of these components. The structure function 
for the magnetic field differs by less than 10 percent from its asymptotic value at a sepa- 
ration of 2/kram- Here, /Cmin is the cutoff wavenumber introduced to prevent an increase in 
the spectrum at wavenumbers that are smaller than those at which energy is believed to be 
injected into the gas (Wallin et al. 1998). We thus adopt the convenient expression 2/kuan 
as an excellent approximation for the correlation length. With the customary assumption 
that alignment of the dust grains is independent of the strength of the magnetic field, the 
behavior of the linear polarization depends on the ratios of the strengths of the random to 
the uniform magnetic fields 

h = -Brms/ -Bavg, (2) 

where -Bavg uniform field component in ith direction, and i represents the field com- 

ponent either perpendicular to the line of sight (_L) and hence in the plane of the sky, or 
it represents the field component that is parallel to the line of sight (||). When there is no 
uniform field in a given direction, hi — oo. Note that the average of the turbulent component 
is zero. 

The fractional linear polarization of an emergent ray of radiation can be expressed as 

when the optical depth and fractional polarization both are small (e.g., Wardle & Konigl 
1990). Here, is a factor that accounts for the imperfect grain ahgnment, is the column 
density of the dust grains, C± and Cy are absorption cross sections perpendicular and parallel 
to the grain symmetry axis, and the effective polarization cross section Cpoi is given by 

^ f Cj_ - C|| (oblate grains) 

~ I - C±) (prolate grains). ^ ^ 
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The relative Stokes parameters are defined as 

q — y cos^ 7' cos 20 ds (5) 

and 

u = y n^cos^7'sin20ds. (6) 

In these expressions, is the number density of ahgned dust grains, 7' is the inchnation of 
the total magnetic field to the plane of the sky at a given point of the ray, and is the angle 
between the projection of the total magnetic field onto the plane of the sky and a reference 
direction — chosen here to be the direction of the projection of the uniform magnetic field 
onto the plane of the sky. If = 0, 7' = const, and does not vary along the ray, eq. (3) 
reduces to 

^-i?^7%rcos^'. (7) 
Lee & Draine (1985) expressed the polarization reduction factor as 

$ = i?FLDCOs2 7, (8) 

where 

to account for the polarization reduction due to a turbulent component of the magnetic field 
(-^ld), due to the imperfect grain alignment (R), and due to the inclination of the uniform 
magnetic field to the plane of the sky (cos^ 7). In eq. (8) 



(picas' 0)-l^ (10) 



where the statistical average is over angles 6 between the local magnetic field and the direction 
of the uniform field. However, the reductions due to turbulence and the inclination of the 
uniform field cannot, in general, be separated as in these expressions. The focus of our work 
is on computing the general polarization reduction factor F due to turbulence, defined here 
as 

F = v/g2 + ^VNd. (11) 

to replace Fld cos^ 7 in equation (8). 

The numerical integration in equation (5) and (6) is performed along the 128^ straight- 
line paths for rays that emerge perpendicular to the surface at the grid points of our "com- 
putational cubes" in which magnetic fields are generated at the 128^ gridpoints. Results for 
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F for the 128^ rays are then averaged to obtain the mean value of F, its standard deviation 
(Tp, and the standard deviation of the polarization position angle cTq, that are appropriate 
when the beam size can be treated as "infinitesimal" (the position angle is obtained from 
the standard relationship a = 0.5tan~^(-u/g). Alternatively, for finite beam sizes, the q and 
u are averaged separately over an appropriate number of rays before F and the direction of 
the linear polarization are calculated. 

3. The Poleirization Reduction Factor F due to Turbulence 

In this Section, we examine the influence of irregularities in the magnetic field on the 
linear polarization of the thermal emission by dust under idealized conditions. Only the 
turbulent magnetic field is allowed to vary along the ray. The density (except for the polar- 
izations based on MHD simulations for the fields) , the degree of alignment and the emission 
by individual grains is taken to be constant and the medium is assumed to be optically 
thin. This model obviously deviates from real clouds which are not uniform and isothermal, 
especially in regions of ongoing star formation. In real objects where there are temperature 
variations, dust emission cannot be presented as emission by a single component (e.g., Vail- 
lancourt 2002). There is evidence that temperature variations in clouds can be modest (e.g., 
Tafalla et al. 2004). Draine & Weingartner (1996) have proposed that radiative torques 
contribute significantly in grain alignment. This suggests that the polarizing power of dust 
grains also varies in real molecular clouds, being favored in regions exposed to the radiation 
field of embedded stars or to the interstellar radiation field. However, in dark molecular 
cloud cores, such as those mentioned in the Introduction, this may not be so important. 
More critical for us are possible changes of the magnetic field parameters that are directly 
or indirectly related to variations in the density. These are neglected in our idealizations. 
Some of these issues will be addressed to a limited degree in Section 4. 

3.1. Infinitesimal Beam Size 

In this subsection, we consider an "infinitesimal" beam size, i.e., the size of the beam is 
much smaller than the characteristic scale length of the turbulence — the correlation length. 
We first treat the case where there is no component of the magnetic field that is parallel to 
the line-of-sight (LOS). The behaviors of the polarization reduction factor F averaged over 
the 128^ rays and the dispersion a„ of the polarization directions of these rays are shown in 
Figure 1 as a function of cloud thickness (in correlation lengths) for clouds with two values of 
b± {b\\ = 0) for which we also have fields from MHD simulations. The polarization is seen to 
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decrease from that due to a completely uniform field (F = 1). As the number of correlation 
lengths through the cloud increases, the polarization approaches the asymptotic value 



of Lee and Drainc which depends only on the ratio of the random to the average magnetic 
field (also noted by Novak et al. 1997). 

The number of correlation lengths needed to reach this asymptotic value varies with bj_. 
If the uniform field is stronger than the turbulent field {bj_ = 0.6, hereafter the strong field 
case), the computed F reaches F^o for a cloud thickness of about ten correlation lengths. For 
the weaker uniform field {b± = 1.5 hereafter the weak field case), iVcorr ~ 20 — 30 is needed 
before the resultant polarization is close to F^. Note that the b± — 0.6 has been designated 
as "medium" in Wiebe & Watson (2001). 

The corresponding dispersion in the position angle is shown in right panel of Figure 1. 
It also depends on the number of correlation lengths that have been traversed and on the 
magnetic field ratio, as shown by Myers & Goodman (1991). However, their expression (the 
dotted curves in Figure 1) gives results that are similar to the computed aa only for uniform 
fields that are rather strong and is less accurate when the uniform magnetic field is somewhat 
weaker than the turbulent magnetic field. 

To justify our procedure for creating the magnetic fields, we also present polarizations 
and dispersions that are computed with magnetic fields which are the result of numerical 
simulations by others for compressible MHD turbulence (Stone, Ostriker, & Gammie 1998). 
They are shown with open and filled triangles in Figure 1. In these simulations, the cor- 
relation lengths and the magnetic field ratios are the same as in the fields created by the 
statistical samplings that are used for the other computations in this Figure. An important 
difference is that in computing F and aa, the variations in the density of the dust along the 
rays are included when we use the results of the MHD simulations. The polarizations ob- 
tained with the weaker uniform magnetic field are in better agreement with those calculated 
with our statistically created fields than are those obtained with the stronger uniform field, 
where larger variations in density along the ray tend to reduce the polarizing power of the 
medium. 

The simphfication that the uniform magnetic field is parallel to the plane of the sky is 
relaxed in Figure 2, and a non-zero line of sight (||-component) is allowed. Specifically, 
we compute F for the same two values of b±^ as in Figure 1 with a range of finite values for fty . 
Note that aa does not depend on b\\ so that the anglular dispersions aa for the computations 
in Figure 2 are given by the right-hand panel of Figure 1 for the same b±. 
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The solid lines with filled and open circles correspond to the two cases shown in Fig- 
ure 1, where the uniform magnetic field along the line-of-sight is zero. When s|vg ^ 0, the 
asymptotic value depends upon fty as would be expected from equation (8), 



However, as can be seen from Figure 3 where asymptotic values arc plotted versus 
for some representative values of cos^ 7, the dependence of the computed F^o on Q and 7 
does not, in general, separate cleanly as in equation (13). In fact, equation (8) only seems to 
express the dependence of F^ on the hne-of-sight magnetic field when < 0.5. For > 1, 
Foo in Figure 3 becomes essentially independent of the hne-of-sight magnetic field. In other 
words, if the turbulent magnetic field is only a factor of 1.5 greater than the component 
of the uniform magnetic field in the plane of the sky, F does not depend on the LOS 
magnetic field component, whether it is weak or strong. This means that our conclusions 
about the polarization properties of the dust thermal emission are essentially insensitive to 
the LOS magnetic field component for the weak field case (6^ = 1-5) that is considered, 
provided 7 < 70°. The number of correlation lengths needed to reach an asymptotic value 
is aproximately the same for any fey. 



To examine the influence of the flnitc size of a telescope beam, we smooth the computed 
map for 6|| = to represent various angular resolutions. As might be expected, the important 
parameter in this case is not the ratio of a beam size to the extent of the map, but the ratio 
S of the beam size to the correlation length Leon- This is illustrated in Figure 4, where 
F and 0^ are presented for three cases in which the ratio of the beam size to Leon is the 
same [S = 1.6), though the dimensions of the computational cubes as measured in terms 
of correlation lengths are different. The areas of the cube surface that arc viewed by the 
beams are adjusted appropriately. Note that Leon is expressed in Figure 4 and elsewhere as 
a fraction of the length of an edge of a computational cube. 

In Figure 5, F and cTq, are presented for several representative beam sizes as measured 
in correlation lengths. When the beam size is much smaller than I/con, F a-nd cTq, are the 
same as for the case of the infinitesimal beam. However, if the beam size is larger than 
the correlation length, the requirement that there be many correlation lengths through the 
cloud in order for F to reach F^o is relaxed. The observed polarization can be close to its 
asymptotic value even if iVcorr < 10. 
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Prom the point of view of obtaining inferences about the strength of the magnetic field, 
this means that low-resolution (or smoothed high-resolution) mapping is preferable. The 
procedure described by Novak et al. (1997) for the large iVcorr limit can then be applied to 
determine the ratio of the irregular-to-average components of the magnetic field even if Neon 
is relatively small or unknown. This straightforward way to infer the irregular-to-average 
magnetic field ratio is valid only if the average magnetic field is not inclined significantly 
to the plane of the sky (that is, 6j| must be small). A large ||-component of the magnetic 
field greatly complicates the interpretation. In fact, an inspection of Figures 2 and 5 shows 
that the curve in Figure 2 for a very strong magnetic field {b± — 0.6 and b\\ — 0.15) looks 
quite similar to the dotted curve in the upper left panel of Figure 5 which corresponds to 
a uniform magnetic field that is weaker than the turbulent magnetic field, but is observed 
with a large beam. 

The rms angle aa decreases with increasing beam size, and even for b± — 1.5 becomes 
smaller than 10 degrees for clouds that are no thicker than about ten correlation lengths 
when (S ~ 3 in Figure 5. This may explain why the measured dispersion of the position 
angle in some observations is comparable to the observational error (Ward- Thompson et al. 
2000). 

For a given value of b± {b\\ — 0), both F and cTa in Figure 5 depend on A^corr and S. 
However, we find that to a reasonable approximation these two quantities can be combined 
into a single quantity Gcorr that describes their efi'ect and can be understood as a generalized 
measure of the number of "correlation cells" that are traversed by rays received by a telescope. 
We define Gcorr as the square of the sum of the square root of the thickness of the cloud and 
the area of the telescope beam expressed in correlation lengths 

G^orr ^ {N'J^r + Sf . (14) 

The validity of this expression is demonstrated in the right-hand-side of Figure 5, where F 
and (Ta are shown as functions of Gcorr- Equation (14) shows that if Ai'corr is comparable to 
S, changes in S tend to be more important for F and aa than are changes in Neon- When 
the number of correlation lengths along a ray is large, Gcorr ~ -Acorr- A rigorous analysis may 
give a more conceptually based relationship for Gcorr- However, the expression (14) serves 
well for the illustrative purpose in this paper. 

Although F and aa each depend on the thickness of the cloud and upon the beam size 
(as measured by S), the plot of F as a function of cTq does not depend upon these two 
quantities. Curves representing F versus a a are shown in Figure 6 for two values of the 
magnetic field ratio b±. We have performed computations with a number of choices for Acorr 
and S to verify that the relationship between F and aa is, in fact, independent of these 
quantities. 



- 11 - 



Also shown in Figure 6 are observational data for a selection of large molecular clouds 

and individual globules (both with and without embedded sources). The data for OMC-1, 
M17, W51, and W3 are taken from Dotson et al. (2000) and Schlcuning et al. (2000). 
Data for NGC 7538 are taken from Momosc et al. (2001). Data for OMC-3 arc taken 
from Matthews et al. (2001). For each object, the average polarization p and the "true" 
dispersions of p and a are computed from 

aj = varp-l/5](l/<jy (15) 

and 

where (Xp^i and a^^i are the observational errors that arc provided by the observers. We exclude 
points with p/ap < 3 as unreliable, as well as points that apparently (based on the position 
angle histogram) belong to a component that is different from the "main" component. Data 
for CB 26, CB 54, and DC 253-1.6 are taken from Henning et al. (2001). As no tables 
are given in the latter paper, all the needed numbers are extracted from the text. Data for 
L1544, L43, and L183 are based on the paper by Crutcher et al. (2004) and kindly provided 
by its authors. 

The wavelengths of the observations are indicated in the caption. It is well known that 
polarization is wavelength-dependent. Hence, the data points obtained at different A cannot, 
in principle, be compared directly, unless they are reduced to a single wavelength with some 
correction factors like those implied by Figure 15 from Hildcbrand ct al. (2000). However, 
there may be considerable spread in these factors from object to object (Vaillancourt 2002), 
so we choose to make no wavelength correction of any kind. In the relevant range of A 
(100 — 850 fim) these factors are near unity. 

To relate the observed percentage polarization to the F that is computed, it is necessary 
to know the polarization that is produced by the grains when the magnetic field is completely 
uniform. Hildebrand & Dragovan (1995) estimated that the maximum polarization produced 
by the mixture of perfectly aligned silicate and graphite grains {R = F = 1) should be about 
35% at A = 100 /im for their particular grain shape. We assume that the reduction from 
the theoretical maximum of 35% to the maximum polarization that is observed (~ 15%) 
occurs entirely because of imperfect grain alignment and there is no contribution due to any 
irregularity of the magnetic field {R ~ 0.4). Thus, to compare the observational data with 
the computed F, we adopt Fobs — Pohs/{R -35%). 

The data points in Figure 6 lie mostly in the region bounded by the two computed curves. 
This tends to indicate that the uniform magnetic field in these objects is comparable with the 
random field or does not exceed it by much. Of course, there is considerable uncertainty in 
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relating the observations to the computations in Figure 6 because of the poor knowledge of R. 
Without information about R, values of p and cTq, can only be used to provide general guidance 
about the ratio of the irregular to the uniform magnetic field. With this limitation in mind, 
Figure 6 demonstrates the key result — that the observed average polarization properties 
of interstellar clouds are consistent with irregularities in the magnetic field and a uniform 
magnetic field that is relatively weak, as we favored in Wiebe & Watson (2001). 

A non-zero component of the uniform magnetic field that is parallel to the line of sight 
does not affect 0"^, but does tend to reduce F and hence shifts the computed curves downward 
in the left hand panel of Figure 6. As long as this parallel component is similar to or less 
than the perpendicular component of the uniform field, its effect (see Figure 2) will be to 
shift the curve for b± = 0.6 downward by no more than about 0.1 in F and to leave the curve 
for = 1.5 essentially unchanged. The data in Figure 6 would still be consistent with the 
relatively weak uniform magnetic field. 

A quantity that is free from the uncertainties in R is the relative polarization dispersion, 
gf/F. In the right hand panel of Figure 6, this quantity is plotted as a function of cTq. The 
relationship cTq = 28.6°(7ir/F, that would be expected for statistically independent q and 
only holds for small a a- For a a > 10°, the value of op/F is only weakly sensitive to b±. It 
also does not depend on the beam size, though this is not shown explicitly. On the other 
hand, it grows somewhat for large h\\ . Thus, a relative polarization error that is greater than 
unity can be indicative of a significant linc-of-sight component of the magnetic field. With 
the few exceptions, almost all the observational points are concentrated near the computed 
values. Note that crp and aa in Figure 6 are not observational errors, but actual dispersions 
caused by the variations of the magnetic field. 

In Figure 7 we show the maps of polarization vectors that have been used to compute 
the average values presented in the foregoing. Maps are shown for two values of Ncoit — the 
number of correlation lengths traversed along each ray. The maps are smoothed with beam 
sizes of 5" = 1.6 and 3.2. They correspond to dotted and short-dashed lines on the leftmost 
panels of Figure 5. While the mean polarization is close to its asymptotic value on these 
maps, there are regions where F exceeds this value almost by a factor of 2. Even though the 
uniform magnetic field is weak, the polarization shows a regular pattern at Ncott ~ 30, when 
the beam size is comparable to the correlation length, and at A^corr ~ 10, when the beam 
size is greater than Lcorr- 

We conclude that the irregularities in the magnetic field used in Wiebe & Watson (2001) 

to understand the polarized absorption of starlight are not in conflict with the observed 
polarization properties of the thermal radiation that is emitted by dust grains, provided 
that the correlation lengths in the regions that have been investigated are somewhat smaller 
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than the sizes of the telescope beams. It is more difficult to understand how a polarization 
hole can appear in such an environment. Along with the regions of high polarization, there 
are some spots in Figure 7 where F is much smaller than the mean value. Inspection of 
the underlying map for the infinitesimal beam size shows that this is due to the higher than 
average magnetic field tangling in the particular regions. If for some reason disturbances in 
the magnetic fields in actual gas clouds are associated with the density enhancements, such 
regions will be observed as polarization holes. This issue is further addressed in the next 
Section. However, in general, it is unlikely that a polarization holes will be due to this cause 
alone. Prom Figure 1, the decrease in the percentage polarization can be caused either by an 
increase in the number of correlation lengths or by a weakening of the mean magnetic field 
(or both) in a dense region where the hole is observed. However, in the weak uniform field 
case and at A^'corr > 10, which is preferable for our proposed interpretation of the polarization 
of starlight, F already is about as small as it can be (especially, when the finite beam size is 
taken into account), leaving little space for further decrease. 

On the other hand, the magnetic field structure may cause polarization holes if for 
some reason the uniform field is stronger on a periphery of a dense object and weakens 
closer to its center. If the decrease of polarization with growing intensity is caused by 
shorter correlation lengths, one would expect that is smaller in dense cores than in the 
surrounding medium. Interesingly, this is what actually is observed in some cores of Barnard 
1 dark cloud (Matthews & Wilson 2002). 

4. Density Effects 

In the previous section we concluded that the unresolved magnetic field structure cannot 
be a common reason for the polarization hole effect. In this Section we investigate whether 
the polarization hole effect can be due to effects that are related to the variations in density. 

We generate the spatial distribution of dust by the same procedure that we used to create 
the magnetic field. A single component of the vector field is created in the Fourier space with 
the same /cmin as is used to generate the magnetic field. This distribution of Fourier waves 
is inverted to find the distribution in coordinate space. The resultant Gaussian quantity 
is interpreted as the logarithm of the density. The density created in this way has a log- 
normal probability distribution function (PDF) that possesses the desired spatial variations 
as characterized by the correlation length. The volume weighted distribution of the quantity 
y = log(nd/rid), where rid and fid are the actual and spatially averaged densities of dust 



-14- 



grams, is 



1 



iy+\f^\) 



(17) 



exp 



2ct2 



The log-normal PDF seems to be relevant in the isothermal case and is reproduced in many 
MHD simulations (Nordlund & Padoan 1999; Vazquez-Semadeni et al. 2000; Ostriker et 
al. 2001). In Figure 8 we compare the PDFs for our simulated densities with that of the 
MHD simulation by Stone et al. (1998). The corresponding values are given in the 
legend. To assess the importance of the spatial structure of the density field, distributions 
with different numbers of correlation lengths (A^con = 5 and 12) through the computational 
cube are considered. The polarization reduction factor F in this Section thus includes the 
reduction of polarization due not only to tangling of the magnetic field, but also due to the 
various additional factors that are studied here. The correlation length is assumed to be 
independent of density. 



In this subsection we examine the consequences of the assumption that dust grains lose 
their ability to polarize light as the density increases. This might be due to poor alignment 
caused by more frequent collisions with the gas molecules or to the growth of icy mantles that 
make the grains more spherical. A similar suggestion was recently studied by Padoan et al. 
(2001). They showed that it is possible to reproduce the observed anticorrelation between P 
and / under the assumption that grains are not aligned at depths beyond a certain value of 
extinction Ay measured from the edges of the clouds. This assumption implies that radiative 
torques are the primary cause for grain alignment. 

Here, we assume that grains are no longer aligned when the density exceeds a certain 
critical density rin^. That is, in equations (5) and (6) we assume that 



First we consider the density distribution that corresponds most closely to the distri- 
bution from the MHD calculations (labeled 'N' in Figure 8). Intensity maps overlaid with 
the polarization vectors and the F vs / scatter diagrams for the 'N' distribution are shown 
in Figures 9 and 10. The density cutoff is set as 1.7 times the median density. The upper 
panels of the Figures correspond to calculations for which the number of correlation lengths 
along the side of the map A^con ~ 5. In the calculations for the lower panels A^con ~ 12. 
The uniform magnetic field is strong {b± — 0.6) in Figure 9, and is weak {b± — 1.5) in 



4.1. A Density Cutoff for Polarization 




rid < nii 
rid > nir 



lim 



lim 
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Figure 10. It is noteworthy that the densities at almost 30% of the grid points exceed the 
specified critical density and thus do not contribute to Q and U . If we assume that the side 
of the computational cube is 1 pc (so that the diameter of the central concentration in upper 
panels is about a few tenths of a parsec and Lcorr ~ 0.1 pc), the optical depth r ^ 10~^ 
for the 'N' distribution corresponds to a median gas density of approximately 10'^ cm"'^ at 
A = 1.3 mm (where a dust opacity Ki.smm = 1 cm^ g~^ is assumed; Ossenkopf & Henning 
1994). To produce the observed decrease in polarization, dust grains should be unaligned in 
all clumps that are at least moderately dense for these calculation with the 'N' distribution. 
This seems to contradict the fact that highly polarized dust emission is sometimes observed 
from regions that are even more dense. 

Even under these rather extreme conditions, the scatter diagram for F vs I for the weak 
field case is far from the tight appearance actually seen in many observations — at least, at 
long wavelengths. In the strong field case, the tightness of the correlation depends upon 
the specific realization at this niim; i- c. in some statistical realizations of the density the 
scatter can be more significant. If num is increased by only by a factor of a few, the correlation 
between the percentage polarization and the intensity completely disappears. Lower values of 
niiin tend to produce tight correlations in any realization, though the percentage polarization 
is quite low. If we set num equal to the median density, the maximum F (at the lowest 
intensity) is only about 0.1 even for the strong field. The observed slope — a decrease by a 
factor of about 3 in percentage polarization as the intensity grows by an order of magnitude — 
is only reproduced by the strong field and low A^corr- 

In other models of HD and MHD turbulence the density distribution can be wider, with 
the range in densities reaching six orders of magnitude (Padoan et al. 1998; Klessen, Heitsch, 
& Mac Low 2000). Hence, we introduce a wider PDF labeled 'W in Figure 8. With the same 
choice of fiducial parameters, the median density is ~ 4 x 10^ cm~^. The high density tail 
of the distribution extends to ~ 10^ cm~^, albeit the number of high density cells is small. 
Thus, only a tiny fraction of the volume contributes to the dust emission. In such a wide 
distribution, the tight correlation between F and / can be reproduced with a much higher 
value of niim- Maps and diagrams for the wide density PDF and ^ SOnmcd are shown in 
Figures 11 and 12. Only 4% of the grid points do not contribute to the polarized emission 
by this criterion. This seems sufficient to obtain (at least in most realizations) a rather tight 
F and / correlation in the strong field case. The weak magnetic field tends to produce more 
scatter, though F at low intensity is quite high in all four cases that we present — reaching 
value of 0.7 for the strong uniform field and A^con = 5. In all cases except the last, F decreases 
significantly with intensity, in agreement with observations. To summarize, increasing riiim 
leads to more scatter; decreasing produces tighter correlations but smaller F. 
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4.2. 



Thermalization Effects 



It is widely assumed that the emission by dust from the star forming regions is optically 
thin. Nevertheless, unresolved optically thick clumps are sometimes mentioned (e.g. Gull 
et al. 1978; Schleuning 1998; Schleuning et al. 2000) along with other possible causes for 
the decrease in the percentage polarization in the brightest regions of dense clouds. Clearly, 
this mechanism is most likely to be effective in star forming regions where high opacities 
and submillimeter luminosities are expected for Class objects (Andre, Ward- Thompson, & 
Barsony 2000). Optical depths of order of 1 are implied both at the far-infrared (Larsson et 
al. 2000; Mookerjea et al. 2000) and in the submillimeter range (Sandell 2000; Sandell & 
Knee 2001) in the dense parts (cores) of at least some star forming regions. 

We will examine, to a limited degree and only in Figure 13, the consequences of the 
thermalization that occurs when the optical depth approaches and exceeds one. For this, 
the previous expressions must be generalized somewhat. The optical depth becomes 



where tq is the normalization factor, and the expressions for relative Stokes parameters are 



u = tqC ^ J cos^ 7' sin 2(f)ds. 

These expressions allow the emission by dust to be optically thick, but still assume that the 
fractional polarization is small. 

To investigate the effect of large optical depth, we utilize an even wider density dis- 
tribution (labeled 'T' in Figure 8) to calculate F for three different values of tq which are 
chosen so that the maximum optical depth along any ray is Tmax = 1, or 10, or 100. There is 
no density cutoff for grain alignment in this case. Results for the second and the third cases 
are shown in Figure 13. At Tmax = 1, no trend of decreasing polarization with increasing 
intensity is evident. It starts to appear at Tmax = 10, where the properties of the density 
field are still realistic. With the adopted fiducial parameters, at A = 1.3 mm the 'T' density 
PDF results in a few rarified cells (n ~ 1 cm~^) and a few very dense cells (with the peak 
gas density n 2 x 10^° cm~^). The median density for the entire cube is 10^ cm~^, and 
median optical depth is ~ 5 x 10~^ (for the far-infrared waveband, all densities are an order 
of magnitude lower due to the higher opacity). 

However, a slope that is similar to what is observed is only achieved when Tmax = 100 
and with a median optical depth ~ 10~^, which probably is too high for the submillimeter 
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q = tqc ^ J cos^ 7' cos 20(is 
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waveband. Thus, we conclude that thermahzation probably is not important for explaining 
the polarization holes observed at longer wavelength, though it may contribute to this effect 
in the far-infrared. 



4.3. Effects of Correlations between the Density and the Magnetic Field 

The procedure that we use to generate the matter densities does not lead to distributions 
for the matter and magnetic fields that are consistent with one another. Specifically, in the 
previous subsections wc used the same ratio of the turbulent to the uniform magnetic field 
b± for the entire cube, despite the variations in density. In reality, we may expect that the 
spatial distributions of density and of the magnetic field are not independent. In Figure 14 
wc plot the standard deviations of the 3D magnetic field along each ray versus the column 
density in arbitrary units for an MHD simulation (weak uniform field case) by Stone et al. 
(1998). Though the lower part of the plot is more densely populated, a trend is evident. For 
clarity, only each 100th point is shown for column densities less than 400. A similar trend is 
observed in the cube computed with the strong uniform field. Thus, in the available MHD 
simulations, the magnetic field does exhibit a stronger variation along the rays with higher 
optical depth — though the range in density is not wide. 

In our model, the ability to reproduce this correlation between the matter density and 
the magnetic field is quite limited. To evalute its importance semi-quantitatively, we calculate 
a model in which b± is proportional to the local density and is normalized to be equal to a 
specific value at the median density, designated as 6med- The 'N' density PDF is utilized at 
-^corr = 12 (these parameters are chosen to reproduce the results of the MHD simulations as 
closely as possible). There is no density cutoff for the alignment of the grains. In Figure 15 
we show F and / for the calculations with bmed = 1-6 and 6med = 0.5. For the two cases that 
are presented, the magnetic field ratios in most grid cells correspond to those of our weak 
and strong fields, respectively. At the low density end of the distribution, the magnetic field 
is strictly uniform; at the high density tail, there is essentially no uniform field. 

The decrease in F with increasing intensity is present in both cases, though only the 
strong field (6med = 0.6) reproduces the observed slope. If we use the same scaling of b± with 
density and the 'W density PDF, the trend is absent both in the weak and the strong field 
cases. It only reappears at b^ed < 0.2. However, this scaling probably is not appropriate for 
the 'W distribution, which is intended to describe the more advanced evolutionary stages 
of the star forming regions. 
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5. Discussion 

The polarization of starlight due to absorption by dust and the polarization of the 
thermal emission from dust provide two ways to study the structure of the magnetic field 
in star forming regions. However, observations indicate that tracing the direction of the 
magnetic field with polarized starlight can be seriously hindered by some factor or factors that 
prevent the grains in the dark interiors of clouds from influencing the observed polarization. 
In Wiebe & Watson (2001) we argued that a relatively strong turbulent magnetic field within 
dark clouds may sufficiently randomize the orientations of the dust that their contributions 
to the polarization of starlight effectively cancel. If a ray of starlight is initially polarized 
by passing through the general interstellar medium with b± ~ 0.5, and then traverses about 
10 correlation lengths within a cloud where b± is about 3 times greater, its polarization is 
almost unchanged in magnitude and in direction as a result of passing through the cloud. 

In this paper we show that the occurence of such irregular magnetic fields in dark clouds 
is not in confiict with the relatively ordered polarization vectors of the dust emission observed 
both in quiescent dense cores and in the regions of ongoing star formation. The maps of 
polarization vectors that are computed using a medium with relatively strong turbulent 
magnetic fields and with large iVcorr show the regular patterns typically seen in the infrared 
and submillimeter observations. The average percentage polarization within these maps is 
2 — 4% (assuming R 0.4), and is close to what is observed. Heitsch et al. (2001) reached 
similar conclusions about the polarized emission for certain specific turbulent MHD models 
for the medium. Chance misalignment of polarization vectors in adjacent points can cause 
a decrease in the percentage polarization — reminiscent of polarization hole — if this region is 
observed with low angular resolution. 

In the context of the polarization hole effect. Figure 1 suggests that the decrease in the 
percentage polarization can be a result of the larger number of correlation lengths that a ray 
traverses in passing through a dense core versus the number it traverses in the surrounding 
medium. For this, Ncott in the surrounding medium must be small because the polarization 
reduction factor F reaches its asymptotic value quite fast. For the moderately strong uniform 
field {b± — 0.6), F is essentially Foo when iVcorr is approximately a few. In the case of the 
weak uniform field {b± — 1.5) that we favored in Wiebe & Watson (2001), a factor of a 
few decrease in F is achieved between A^con ~ 1 and 30. Over the same range of A^con, c^a 
decreases from 40° to about 20°. 

These values for A^con and ctq, are larger than the values that are indicated from the 
observations. The need to have A^con near 30 is relaxed if we take into account a finite beam 
size and allow the correlation length Lcorr to be smaller in the core than on the periphery. 
Let us assume that on the periphery of a dense core, a ray traverses A^con ~ 1 — 3 and that 
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the beam size is smaller than Lcorr- In this case F is about 0.4 — 0.5 (the sohd or long-dashed 
curve on the top left panel of Figure 5). If toward the center of the dense core Leon becomes 
smaller than the beam size, then F ~ 0.1 toward the center (the dotted curve on the top 
left panel of Figure 5). With this beam size and with Ncott ~ 10, cTq, becomes less than 10°. 
If the same region is observed with a smaller beam (e.g., with an interferometer), then the 
dispersion in position angles would increase (because Gcorr becomes smaller). 

The apparent discrepancy in this scenario of a large dispersion in position angles at the 
periphery of the core (~ 40°) can be partially circumvented by assuming that the strength of 
the uniform magnetic field is higher in the vicinity of the core than in the core itself. When 
b± = 0.6 in the periphery region, aa can be as small as 20° (solid curve on the lower right 
panel of Figure 5 at Neon = 3). It is even smaller at lower b±. 

The requirement that a ray traverse only a few correlation lengths in a region surround- 
ing the dense core where the polarization hole is observed may also seem to conflict with our 
conclusion in Wiebe & Watson (2001) that a typical dark cloud spans no fewer than 10 cor- 
relation lengths. It must be noted that in the Bl cloud, where polarimetry observations of 
both starlight and thermal dust emission are available, they indicate different directions for 
the magnetic field (Matthews & Wilson 2002). It is quite possible that the polarized hght of 
background stars and the polarized dust emission, even when they are observed close to one 
another, trace different objects or parts of the same object that differ in physical conditions 
and/or scale. The two types of observations may thus not be directly related to one another. 
In addition, observations of starlight polarization arc mostly used to study quiescent dark 
clouds, while many (though not all) objects in which the polarization of the dust emission 
is observed already contain Class protostars with outfiows. 

Our calculations demonstrate quantitatively how the magnitude of the correlation length 
(or at least how the correlation length compares with the size of a cloud and the size of the 
telescope beam) is the key factor for the inffuence of turbulence on the observations of the 
linear polarization. Unfortunately, any suggestions based on our calculations that we can 
offer for extracting the correlation length from observational data are hampered by the lack 
of information about R and about the source geometry. With these limitations in mind, 
some inferences about iVcorr can be made based on Figures 5 and 6 and an application of 
equation (14). Assuming that the number of correlation lengths is the same along the line 
of sight and perpendicular to it, we may replace 5* in equation (14) with Ncorr/F), where D is 
the number of beams across an object. Equation (14) then can be solved for Neon in terms 

of (j^corr* 

Let us consider two examples. The location of the L1544 core in left panel of Figure 6 
implies that the uniform magnetic field is strong in this object. The polarization position 
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angle dispersion ctq then corresponds to Gcon ~ 10 (see right panel of Figure 5). For these 
observations, D ^ 8 and we obtain Ncorr ~ 5. Assuming the generally accepted distance of 
140 pc and an angular size of about 110" (Ward-Thompson et al. 2000) gives L^orr ~ 0.02 pc 
for this object. The L43 core, on the other hand, in left panel of Figure 6 is located close 
to the curve corresponding to the weak uniform field. The small anglular dispersion implies 
Gcorr ~ 50 according to Figure 5 (right panel), which for the same D gives iVcorr ~ 20. 
At a distance of 170 pc and an angular size of ~ 100" (Ward-Thompson et al. 2000), 
-^corr ~ 0.004 pc. Determinations of the correlation lengths by independent methods would 
be valuable. Unfortunately, the tendency for velocity shifts in clouds to be caused by large- 
scale velocity gradients as well as by turbulence, for the presence of inhomogeneities in cloud 
properties, etc. obscure the detailed effects of turbulence . 

This work was supported in part by NSF Grant AST 9988104. DW acknowledges 
support from the President of the RF grant MK-348.2003.02. 
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Fig. 1. — Polarization reduction factor F due to turbulence and the dispersion in the 
position angle vs the number of correlation lengths A'^con through a cloud. Results from the 
actual MHD computations are denoted by filled {b±^ = 1.5) and open {b± = 0.6) triangles. 
The dispersion of position angles obtained with the expression of Myers & Goodman (1991) 
is indicated by the dotted lines. 
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Fig. 2. — Polarization reduction factor F due to turbulence and with a non-zero LOS mag- 
netic field. Weak and strong field cases are indicated by open and filled circles, respectively. 
Solid lines = oo) are computed when the uniform magnetic field component along the 
LOS is zero. 




Fig. 4. — Polarization reduction factor F due to turbulence and the dispersion cr^ in the 
position angles vs the number of correlation lengths iVcorr through the cloud for various values 
of the correlation length and h± = 1.5. The ratio of the beam size to Lcorr is the same for 
all three sets of curves. 
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Fig. 5. — Left panel — F and cTq, vs the number of correlation lengths Ai'corr through a cloud 
for different magnetic field ratios and beam sizes. Right panel — Same quantities as in left 
panel, but shown as functions of a generalized number of correlation lengths Gcorr- 
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Fig. 6. — Computed polarization reduction factor F due to turbulence (left panel) and the 
relative polarization error ap/F (right panel) as a function of cTq,. Observational data also 
are shown for comparison. The name of the object is centered on the corresponding data 
point. The number 100 or 350 next to the name of the object indicates the wavelength in 
fim. If no number is specified, the measurement is at 850 fim. Some cloud labels have been 
shifted slightly for clarity 
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Fig. 7. — Computed maps of polarization vectors with the weak, uniform magnetic field for 
two values of N^orr (left column, Ncorr = 12; right column. Neon = 36) and for two values of S 
(top row, S = 3.2; bottom row, S = 1.6). Coordinates are given as a fraction of the distance 
along the face of the computational cubes. The longest vector corresponds to F = 0.18. The 
uniform field is in the x direction. Position angle dispersions are 8° (left) and 7° (right) for 
the top row, and 21° (left) and 11° (right) for the bottom row. 
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Fig. 8. — Density distributions that are utilized in the computations. The distribution 'N' 
('narrow') is chosen to be as close as possible to the results of Stone et al. (1998), which are 
indicated by the thick gray line. Distribution 'W ('wide') is intended to represent results 
of other numerical simulations that lead to a wider density PDF. The distribution used to 
study thermahzation effects is labelled with 'T'. 
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Fig. 9. — Map of polarization vectors and F ys I diagram for the 'N' distribution and the 
strong uniform magnetic field. Upper panel, iVcorr = 5, lower panel, iVcorr = 12. The longest 
vector corresponds to the highest F. The map is obtained by averaging over 16 x 16 emerging 
rays. The grayscale on left panel and the scale of the horizontal axis on right panel is the 
intensity of dust emission in arbitrary units. The uniform magnetic field is in the x direction. 
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Fig. 10. — Same as in Figure 9, but for the weak uniform magnetic field. 
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Fig. 11. — Map of polarization vectors and F ys I diagram for the 'W distribution and the 
strong uniform magnetic field. Upper panel, Ai'corr = 5, lower panel, Ai'corr = 12. Intensity is 
in arbitrary units. The longest vector corresponds to the highest F. The map is obtained 
by averaging over 16 x 16 emerging rays. The grayscale on left panel and the scale of the 
horizontal axis on right panel is the intensity of dust emission in arbitrary units. The uniform 
magnetic field is in the x direction. 
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Fig. 12. — Same as in Figure 11, but for the weak uniform magnetic field. 
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Fig. 13. — Polarization reduction due to effects of thermalization. 
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Fig. 14. — The rms magnetic field versus column density in the MHD simulations by Stone 
et al. (1998), in arbitrary units 
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Fig. 15. — Polarization reduction factor F as a function of intensity for variable 6j 



